Clustering in the nuclear Fermi liquid 
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We analyze the emergence of various structures in nucleonic matter, such as crystal, clusters, 
liquid drops and haloes. The formation of clusters indicates that nuclei behave like a Fermi liquid 
close to the liquid to solid transition. The relevant parameter is the ratio of the dispersion of the 
single-nucleon wave functions in the nucleus to the inter-nucleon distance. We also discuss the 
relationship between cluster states in nuclei and the pasta phase in the crust of neutron stars, as a 
transitional state between a Fermi liquid and a crystal. Haloes and clusters exhibit opposite features 
with respect to nucleonic localization. 
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Nucleons in atomic nuclei and extended nuclear mat- 
ter exhibit a variety of phases. Liquid drop aspects, for 
instance, were first inferred [lj from fission properties in 
heavy nuclei. Soon after it was also predicted and ob- 
served that cluster states could occur, especially in light 
nuclei Q- Halo structures in nuclei were discovered in 
the late 1980's Q. Although various theoretical mod- 
els have been developed that successfully describe the 
corresponding characteristics of theses states, open ques- 
tions remain: can nucleonic crystal states occur?; do all 
nucleonic phases (liquid, cluster, halo, crystal) have a 
common origin and, therefore, can they be described in 
a unified theoretical framework? In a recent study [J] we 
have shown that the confining nuclear potential deter- 
mines the degree of clustering in nuclei. In the present 
work we aim to further provide a unified view of the var- 
ious nuclear states: cluster states should be considered 
as a transitional phase between a quantum liquid (nu- 
clear matter) and a solid (crystal) and, in general, nuclei 
can be considered as liquid drops close to crystallization. 
These considerations are also relevant for the description 
of the crust of neutron stars, where it is known that de- 
creasing matter density (further from the center of the 
star) leads to a transition from the nuclear matter phase 
(liquid) to a Wigner crystal, with a pasta (cluster) phase 
in between 0-0| . 

Clustering - the arrangement of nucleons in clusters 
of bosonic characters, especially in light nuclei, coexists 
with the nuclear mean- field. The mechanism that drives 
the transition from the fermionic liquid to the crystal-like 
cluster structures is not fully understood [lj but, as we 
have shown in the recent study of Ref. [3] , its origin partly 
lies in the nuclear potential that confines the neutrons 
and protons to the nucleus. The deformation of a nucleus 
also plays an important role in the formation of clusters 
because it removes the degeneracy of single-nucleon lev- 
els associated with spherical symmetry [2j. Clustering 
predominantly appears in excited states, and this can be 
understood from the fact that closeness to the particle 
emission threshold favors cluster formation. States close 
to the continuum cannot be isolated from the environ- 



ment of scattering states, so cluster states at the thresh- 
old belong to an open quantum system The origin 
of cluster formation, however, lies in the effective nuclear 
interaction, and a fully microscopic description of cluster- 
ing necessitates a framework that encompasses both clus- 
ter and quantum liquid-drop aspects in light and heavier 
nuclei. In Ref. [|| we performed self-consistent mean-field 
calculations based on nuclear energy density functionals 
(EDFs) to explore the effect of the confining nuclear po- 
tential on the formation of localized structures in the 
equilibrium state of 20 Ne, and showed that the depth of 
the potential drives the energy spacings between single- 
nucleon orbitals in deformed nuclei, the localization of 
the corresponding wave functions and, therefore, the de- 
gree of nucleonic density clustering. The localization of 
the single-nucleon wave function can be quantified by the 
dimensionless parameter Q: 
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where f is the average inter-nucleon distance and Ar the 
spatial dispersion of the wave function. When a > 1 nu- 
cleons are delocalized and the nucleus displays a quantum 
liquid structure. On the other hand, when a < 1 nucle- 
ons are localized on the nodes of a crystal-like structure. 
The transition from the quantum liquid phase to cluster 
states occurs for a ~ 1. The s, p and d wave functions 
give the principal contribution to cluster states in light 
nuclei. The spatial dispersion of can be calculated from: 
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In Ref. [4( we used two functionals that are representa- 
tive of the two standard classes of nuclear EDFs: the non- 
relativistic Skyrme SLy4 [l(|, and the relativistic func- 
tional DD-ME2 [11], to calculate the self-consistent equi- 
librium mean-field solution for 20 Ne. Even though both 
functionals reproduce the empirical ground-state proper- 
ties (binding energy, charge radius, matter radius) with a 
typical accuracy of 1%, as well as the quadrupole defor- 
mation of the equilibrium shape, the density calculated 
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Self-cons. 


HO 


SLy4 


DDME2 


SLy4 


DDME2 


20 Ne 


0.987 


0.971 


0.996 


0.973 


24 Mg 


1.001 


0.948 


1.015 


0.963 


28 Si 


0.986 


0.955 


1.053 


0.996 


32 g 


0.991 


0.961 


1.057 


0.993 


208 p b 


1.278 


1.309 


1.458 


1.395 



TABLE I. Left: the localization parameter a calculated 
from self-consistent equilibrium solutions determined using 
the SLy4 and DDME2 (TJ density functional (Eqs. Q 
and ([2]), see text). Right: the same but using the 3D harmonic 
oscillator approximation (Eq. (B}). 



with SLy4 displays a smooth behavior characteristic of 
a Fermi liquid, whereas the functional DD-ME2 predicts 
an equilibrium density that is much more localized, with 
pronounced cluster structures. The difference in the cal- 
culated equilibrium density distributions was traced back 
to the depth of the confining potential: in the centre 
of the nucleus the depth of the DD-ME2 single-neutron 
potential is -79.7 MeV, whereas the depth of the SLy4 
potential is -71.4 MeV. The corresponding values of the 
single-proton potentials are: -74.8 MeV for DD-ME2, and 
-66.5 MeV for SLy4. 

Using the same functionals SLy4 and DD-ME2 to cal- 
culate the self-consistent equilibrium nuclear mean-field 
state, in the present work we have computed the spatial 
dispersion of the Nilsson state [110 1/2] in the N = Z 
nuclei 20 Ne, 24 Mg, 28 Si, 32 S, using Eq. ©. The disper- 
sion of the same state is also calculated for 208 Pb. Taking 
f = 0.9 fm as a characteristic inter-nucleon equilibrium 
distance [lH, we determine the corresponding values of 
the localization parameter a displayed in the first two 
columns of Table |U In light nuclei the values of a calcu- 
lated with DD-ME2 are smaller than those obtained us- 



ing SLy4, reflecting the more pronounced localization of 
the nucleonic densities that was already observed in our 
previous study in Ref . [3| • While for light nuclei a < 1, in 
the case of 208 Pb a is considerably larger than 1 and this 
unambiguously characterizes the quantum liquid phase 
of nucleonic matter in this nucleus. 

In a non self-consistent approximation one can use a 
3-dimensional isotropic harmonic oscillator (HO) for the 
the confining nuclear potential. This approximation al- 
lows for a further qualitative discussion of the effects of 
the effective nuclear interaction on the spatial arrange- 
ment of nucleons. The 3-dimensional HO wave functions 
tpkimif) for the first s, p and d states read [13j: 



where b is the oscillator length defined by: 
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muj o (2mVb) 1/4 ' 
R is the radius of the potential for which V = 0, and 
Vo denotes the depth of the potential for r — 0. We 
note that charge radii of atomic nuclei are determined 
with high precision in electron scattering experiments, in 
contrast to the depth of the confining potential Vb which 
is experimentally poorly constrained. 

A straightforward calculation using Eqs. © and Q 
yields Ar 0.56 for the first s, p and d states, and they 
display a Gaussian-like radial dependence (Eq. ([3])). Con- 
sequently, for a constant radius R a deeper potential Vo 
implies a smaller value of the oscillator length b (Eq. (j4])), 
and thus a smaller dispersion. This effect is illustrated 
in Figure 3 of Ref. [H , which shows the radial part of the 
p wave function for two different values of the oscillator 
length b. A smaller value of b implies a more localized 
wave function, eventually leading to cluster formation. 

This concept can be extended to the more general case 
of deformed nuclei by approximating the confining po- 
tential with an axially deformed HO. The wave functions 
are then expressed as [3, EH 
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where H and L are the Hermite and Laguerre polyno- 
mials, respectively. Eq. ([5]) displays a radial dependence 
similar to the 3D isotropic case (Eq. the dispersion 
of the wave functions depends on the oscillator lengths 
b z and bj_ in the respective directions, which in turn de- 
pend on the depth of the potential. The effect of the 
depth of the potential on the localization of the nucleon 
wave functions is universal. 

In the second two columns of Table Q] we list the val- 



ues of a obtained using expression (f?]) for the harmonic 
oscillator length: 



r r (2mV ) 1/4 ' 

with R and Vq determined microscopically using the SLy4 
and DD-ME2 self-consistent equilibrium solutions, re- 
spectively, and ro = 1.25 fm [l6j]. The same trend is 
found as in the self-consistent calculation: a is smaller 
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FIG. 1. The localization parameter a as a function of the 
number of nucleons (Eq. [7])). The transition to cluster states 
occurs for a ~ 1. 



in the case of DDME2, and a > 1 for 208 Pb, indicating 
a quantum liquid behavior for this heavy nucleus. 

The localization parameter a generally increases with 
the number of nucleons A: cluster states are more likely 
formed in light nuclei, as observed experimentally [2j. 
Using the liquid drop parameterization for the radius 
R=r A 1 / 3 , Eq. (JBJ reads 
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Figure [T] displays the evolution of the parameter a with 
A, for typical values of Vo= 70 MeV. Obviously clus- 
tering is favored in light nuclei, whereas the quantum 
liquid phase is predominant in heavy nuclei. Quantita- 
tively, the transition from a cluster to a liquid state in 
nuclei is obtained around A « 30 (a between 0.8 and 
1.2), and this is in agreement with the empirical obser- 
vation that cluster states are more likely to occur in light 
nuclei. For heavier systems a is considerably larger than 
1 and, therefore, cluster states are less likely to appear 
in heavy nuclei. It should also be noted that Fig Q] nicely 
illustrates why a crystal phase (a at least below 0.8) can- 
not occur in nuclei. However, as described below, nature 
may offer the possibility of existence of nucleonic crys- 
tals in neutron stars, due to gravitational and density 
conditions at work at these astrophysical sites. 

The issue of solid (crystal) vs. quantum liquid nature 
of nuclei was already addressed by B. Mottelson, who 
introduced the quantality parameter : 



A= 



mr 2 Vq 



(8) 



where V'o is the strength of the bare nucleon-nucleon 
interaction (~ 100 MeV), and f is inter- nucleon equilib- 
rium distance. The quantality A is defined as the ratio 



of the zero-point kinetic energy to the potential energy, 
and typical values obtained for nuclei are characteristic 
for a quantum liquid phase [Hj. However, cluster states 
in nuclei introduce an additional phase of nucleonic mat- 
ter. The parameter a introduced in this work appears 
more accurate in this respect than the quantality be- 
cause a takes into account the nuclear medium effect on 
the potential, whereas A uses the one that corresponds 
the bare nucleon-nucleon interaction. The quantality 
may not contain information of the localization of the 
wave-function in the nuclear medium. Two characteris- 
tics lengths (Ar and f) should be considered to qua ntify 
crystalline vs. liquid effects, as in condensed matter (171 ]: 
the localization of the constituent wave function in the 
system, and the inter-constituent distance. 

An important ingredient to understand clustering is 
the repulsive part of the interaction which gives rise to 
a crystal lattice on which nucleons can be localized. For 
instance, in the crust of compact stars the Coulomb force 
plays the role of the repulsive interaction that frustrates 
the system and leads to nucleonic crystallization in the 
external low-density part of the crust. The so-called 
pasta phase, predicted in the inner part of the crust is 
considered as a transition between a liquid and a Wigner 
crystal phase [5|, |6[ . An analogy can be drawn between 
nuclear cluster states and the pasta phase for which the 
molecular dynamics approach is used |7j . Interestingly, in 
the seminal paper on pasta phase Baym et al. discussed 
about the necessity to investigate the transition between 
the solid and liquid phases in neutron star matter 
Lattimer and Swesty considered the ratio u of the vol- 
ume of a nucleus to the volume of the Wigner- Seitz cell 
Q. This quantity corresponds to the cube of the ratio 
of the nuclear size to the inter-nuclear distance, and this 
is analogous to the definition of a (Eq. ([1])). In Ref. 
6] u — 1 corresponds to nuclear matter and the pasta 
phase can already appear for u = 0.72. For lower values 
of u the crystal structure dominates. The transition be- 
tween the Wigner crystal and the quantum liquid in the 
neutron star crust can be described by various models: 
gelification [TJ|], Coulombic frustration [ll| or quantum 
melting [2fjj . 

In the case of atomic nuclei the inter-nucleon interac- 
tion contains both an attractive and a repulsive com- 
ponent. Therefore the hard-core part of the nucleon- 
nucleon interaction generates a typical f spatial scale. 
Most nuclei, of course, have a > 1 in their Fermi liquid 
states, especially heavy ones as discussed above. Nuclei, 
therefore, consist of largely delocalized nucleons and this 
explain their liquid drop nature and also the large mean 
free path of nucleons. But, more precisely, nuclei are 
in the Fermi liquid state at the limit of crystallization, 
explaining why localized cluster states can be formed in 
light nuclei. The short range of the hard-core also ex- 
plains why genuine nucleonic crystals cannot occur in 
nuclei, whereas they can appear in the crust of neutron 
stars where crystallization is caused be the long range 
Coulomb interaction. 



4 



Halo states are known to have delocalized wave func- 
tions [2l| and could be related to the present discussion. 
When the confining nuclear potential is approximated 
by a square-well, the oscillations of the wave function of 
a state of energy E (<0) are determined by the wave 
number h 2 k 2 = 2m(E + Vq), whereas outside of the po- 
tential the decay of the wave function is governed by 
h 2 k' 2 — —2mE. Inside the potential a larger k favors 
localization. This happens when the potential is deeper 
and/or E gets closer to 0: the degree of localization de- 
pends on the difference E — (—Vq), and a deeper poten- 
tial favors the localization of the wave function. Also 
a weakly-bound state will be more localized, and this 
means that cluster states are more likely to occur for 
weakly-bound states close to the emission threshold, in 
agreement with the Ikeda picture [22|. An important 
characteristic of nuclei is that quantitatively the disper- 
sion of the single-nucleon wave function can be of the 



same order of magnitude as the inter-nucleon distance, 
leading to clustering as discussed above. Outside of the 
potential the decay of a single-nucleon wave function is 
only determined by the energy of the state, and a weakly- 
bound state displays a large radial extension in the ex- 
ternal region. In the limiting case of very weakly bound 
states, the pronounced derealization of the correspond- 
ing wave functions outside of the potential gives rise to 
halo structures [2l[ . Inside the nuclear well, as discussed 
above, both the energy of a state and the depth of the po- 
tential play a role in the localization of the wave function. 
To summarize, the single-nucleon equation of motion re- 
lates energies (depth of the potential and/or binding en- 
ergies) to the localization of the corresponding wave func- 
tion. In nuclei the typical values of length scales and en- 
ergies allow the formation of clusters, liquid drops and 
haloes. 
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